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We discuss the properties of an accretion disk around a star with parameters typical of 
classical T Tauri stars (CTTS), and with the average accretion rate for these disks. The disk is 
^ I assumed steady and geometrically thin. The turbulent viscosity coefficient is expressed using 

the a prescription and the main heating mechanisms considered are viscous dissipation and 
irradiation by the central star. The energy is transported by radiation, turbulent conduction 
and convection. We find that irradiation from the central star is the main heating agent of the 
^ ' disk, except in the innermost regions, R < 2 AU. The irradiation increases the temperature of 

, the outer disk relative to the purely viscous case. As a consequence, the outer disk (i? > 5 AU) 

' becomes less dense, optically thin and almost vertically isothermal, with a temperature 

\^ • distribution T cx R~^^^. The decrease in surface density at the outer disk, decreases the disk 

, mass by a factor of 4 respect to a purely viscous case. In addition, irradiation tends to make the 

' outer disk regions stable against gravitational instabilities. 
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Introduction 



The co-planarity and circularity of planetary orbits in our Solar System support the idea that it has 
■ been formed from a rotating "protoplanetary disk" , where energy is dissipated and momentum redistributed 

before the time of planet formation. Circumstellar disks are naturally formed during the collapse of a cloud 
fragment with non-zero angular momentum ( passen fc Moosman 1981 ; Tereby, Shu & Cassen 1984). Disks 



produce observational signatures in the spectral energy distributions (SED) of young stellar objects (YSO) 
in the form of excesses of emission at UV, IR and mm wavelengths. 

During the early evolutive phase of a forming star, the disk-star system is still surrounded by an 
infalling envelope, which dominates the far IR SED. At longer wavelengths the infalling envelope becomes 



optically thin, and the SED is dominated by emission from the disk (Kenyon, Calvet & Hartmann 1993 



Calvet et al. 1994; Hartmann, Calvet fc Boss 1996). Radiation from the infalling envelope is one of the 
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most important heating mechanisms of the disk in this early phase and has to be included to understand 



the long wavelength SEDs of embedded stars (e.g. Butner et al. 1991, 1994; DCH97) 



There can be a transition phase in which the star and the disk are surrounded by a tenuous dusty 



envelope which scatters stellar radiation onto the disk (Natta 1993), heating its outer regions. But finally, in 



a more evolved phase, it is expected that the material surrounding the disk-star system becomes negligible 



e.g. Shu, Adams fc Lizano 1987 ) 



Models for accretion disks (Lyndcn-Bell & Pringle 1974) or fiattened disks reprocessing stellar radiation 



(Adams et al. 1987) predict steeper SEDs than those observed in Classical T Tauri stars. Kenyon & 
Hartmann (1987, hereafter KH) have proposed that a disk in vertical hydrostatic equilibrium, with well 
mixed gas and dust, is flared and intercept more stellar radiation than a flat disk. This model predicts a 
SED more similar to the observed ones, but it is restricted to wavelengths A < 100 /im because of their 
assumption that the disk is optically thick. 

In this paper, we explore the evolved phase in which the disk-star system is surrounded by a negligible 
amount of dust, adopting the basic idea behind the KH model, i.e. gas and dust in the disk are well mixed 
and thermally coupled. The disk receives radiation directly from the central star and it is also heated by 
viscous dissipation, cosmic rays and radioactive decay. We calculate the detailed vertical structure of such 
a disk taking into account several transport mechanisms: turbulent flux, radiation and convection. The 
disk interior and atmosphere are not arbitrarily separated. The complete structure, from the midplane to 
the surface, is calculated with the same set of equations, written in a general way to treat both optically 
thick and thin regions. We then present a detailed study of one case, with parameters typical of Classical T 
Tauri stars. 

In a forthcoming paper (paper II) we show results of the structure and physical properties for a 
wide range of parameters characterizing disks and central stars. Finally, motivated by the large amount 
of observational information compiled in recent years, we calculate (paper III) how the observational 
signatures of disks depend on their physical properties. Comparing models with observations of Classical T 
Tauri stars, we infer disks physical properties. 

The structure of the present paper is as follows: in §2 we present the assumptions and general 
description of the model; in §3 and §4 we give a detailed discussion of the heating of the disk and the 
energy transport; in §5 and §6 the equations and boundary conditions are written; in §7 a particular 
case is discussed, with parameters typical of Classical T Tauri disks; in §8 we compare the model with a 
non-irradiated accretion disk and discuss the effects of irradiation; in §9 the effect of the irradiation from 
an accretion ring is calculated and finally §10 presents a summary of the results. 



2. General description of the Model 



We assume that the disks are in steady state, thus, all their properties are time independent. We also 
assume that the disk is geometrically thin, so the radial energy transport is neglected and the vertical and 
radial structures are treated as decoupled problems (Lynden-Bell & Pringle 1974; Pringle 1981). The disk 
is assumed in vertical hydrostatic equilibrium. Since there can be turbulent and convcctivc motions in the 
vertical direction, the assumption of hydrostatic equilibrium refers to the mean bulk motion of the gas and 
dust in the disk. These assumptions are basic for the model and cannot be relaxed without changing the 
disk structure equations and the numerical method used to integrate them. Other additional assumptions, 
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required to quantify the disk physical properties are that the mass accretion rate M is uniform through the 
disk and that the turbulent viscosity is described by the a viscosity prescription. Therefore, the viscosity 
effective coefficient is given by j/j = a c^i?, where Cg is the local sound speed, H is the local scale height of 
the gas and a is a free parameter, which is assumed constant through the disk with the only constrain that 
a < 1 ( 3hakura fc Sunyaev 197^ ). This prescription could be relaxed if M{R) and i't{R,z) are given. 



We consider that the disk is heated by viscous dissipation, radioactive decay, cosmic rays, and stellar 
irradiation. The first two mechanisms generate and deposit energy at each height. Cosmic rays and stellar 
irradiation penetrate the disk from the upper surface, and interact with the material below. We assume the 
disk is a mixture of gas and dust and both components are thermally coupled and described by a unique 
temperature. The gas is locally heated or cooled by viscous dissipation, ionization by energetic particles, 
collisions with dust, dust and gas radiation, convection, and turbulent conduction. The dust grains are 
heated or cooled by radiation from the star and from the disk itself, and also by collisions with gas particles. 
A more general description of the disk structure is required to calculate the disk vertical structure equations 
assuming that gas and dust have different temperatures. However, we assume that the coUisional and/or 
radiative heating between gas and dust are efficient enough at every point in the disk and use only one 
mean temperature to describe both components. 

Because the disk is steady, the energy absorbed is emitted and transported outwards. The energy is 
transported by: (a) a turbulent flux, computed self-consistently with the viscosity coefficient used to describe 
the viscous energy dissipation; (b) radiation, through the first two moments of the transfer equation ; and 
(c) convection, described by the mixing length theory, taking into account that the convective elements lose 
energy by radiation and turbulent flux. 

Given a central star, the important quantities in describing a disk model are M , which gives the total 
energy flux generated by viscous dissipation at each annulus in the disk and a, which quantifies the fraction 
of viscous flux produced at each height. The disk structure also depends on the central star properties, i.e., 
its mass M*, radius i?* and effective temperature T,. The mass is related to the depth of the stellar 
gravitational potential, which is the last source of energy in the accretion process. The stellar luminosity, 
related to and i?*, controls the amount of energy introduced in the disk by stellar irradiation. 

The following sections give the physical and mathematical description of these processes. The resultant 
set of equations are solved with appropriate boundary conditions to get the temperature and density at 
each radius and height of the disk. The disk geometry is shown schematically in Figure |l]. 



3. Heating sources 

The disk heating sources are viscous dissipation, ionization by energetic particles and stellar radiation. 
In this section we describe in detail each of these mechanisms and give the corresponding rates to be used 
in the calculation of the vertical structure. 



3.1. Viscous dissipation 



The energy rate per unit volume locally generated by viscous stress is given by (Pringle 1981 
King fc Rainc 1992|) : 



Frank. 
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(1) 



where p is the mass density, P is the gas pressure, fl is the Keplerian angular velocity, and vt is the 
turbulent viscosity coefficient, given hy vt — a H Cs- 

The total energy per unit area produced by viscous dissipation in a given annulus at a distance R from 
the central star is thus, 



dz = R 



1/2. 



(2) 



where Zoa is the height of the disk, and G is the gravitational constant. 



3.2. Ionization by energetic particles 

The main source of ionization by energetic particles are cosmic rays. However, Stepinski (1992) has 
shown that ionization due to energetic particles produced during the radioactive decay of ^^Al can be more 
important in those regions where cosmic rays can not penetrate. We have thus considered both mechanisms 
and the resulting heating rate per unit volume is written as. 



r,o„(S) = UH, [AQeo. Ccos e-^/^ + AQrad CradC'^ Al)] , (3) 

where 71^2 is the number density of Hydrogen molecules, and A = 96 g/cm^ is the attenuation surface density 
scale of cosmic rays. The ionization rates due to cosmic rays and radioactive decay are Ccos = 10^^^ s~^ 



( iNakano fc Umebayashi 1986| ) and CradC^^Al) = 5.2 x 10"^^ s~i ([Stepinski 1992| ), respectively, and E is the 
vertical column density from the disk surface to a height z. We have taken the energy liberated in heating 
in each cosmic ray ionization as AQcos ~ 20 eV, and for radioactive decay ionization, AQrad ~ lOeV 



( poldsmith fc Langer 197^ ). 



Integrating equation (H) in z, the total energy input per unit area due to ionization by energetic 
particles is given by: 

F^on{Zoo) = / T.on dz - -^{AQ^os Ccos A [1 - g-^-/^^] + AQrad Crad ^} , (4) 

Jo ^ruH 2 

where niH is the mass of the hydrogen atom and Sqo is the total column density of the disk, 

Soo = p dz . (5) 



3.3. Stellar irradiation 

Following Calvet et al. (1991, hereafter CPMD), it is assumed that the stellar radiation intercepted by 
the disk consists of a parallel beam carrying a flux Firr per unit area of the disk surface, incident at an 
angle — cos^^ /io from the normal to the boundary of the disk atmosphere. A fraction of this energy 
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is scattered creating a diffuse field with wavelengths around the characteristic wavelength of the stellar 
radiation, given by the stellar effective temperature T,. The remaining fraction of the incident radiation 
field is truly absorbed and remitted at wavelengths determined by the local temperature of the disk. As 
long as the characteristic temperature of the disk is lower than the effective temperature of the star, it can 
be assumed that no true emission of the disk occurs in the "stellar" frequency range and no scattering of 
the incident beam occurs in the "disk" frequency range. 

The true absorption of radiation at the stellar frequency range is described through a Planck-type 
mean opacity, calculated evaluating the true absorption coefficient at the local temperature and pressure 
but using the Planck function evaluated at T, as the weighting function. The resulting mean absorption 
coefficient is given by: 



/•OO 

K*p{T,P,n)^^^J^ K,{T,P)B,{T^)d,y , (6) 

where is the monochromatic true absorption coefficient, evaluated at the gas temperature T and gas 
pressure P, and is the Stefan-Boltzmann constant. The mean scattering coefficient ap{T, P,T^,), is 
calculated with the same kind of average given by equation (^, substituting k^, by a^, the monochromatic 
scattering coefficient. The mean extinction coefficient is XpiT, P,T^) = Kp + cTp. Since we assume that 
gas and dust are well mixed in the disk at every height, the absorption and scattering coefficients include 
the contribution of dust, unless it is sublimated in regions with a temperature higher than the sublimation 
temperature. If the dust has settled down in the midplane, the upper disk atmosphere would have a lower 
opacity to the stellar radiation than what we are assuming here. 

The impinging stellar radiation field, which propagates in the direction defined by /xo, reaches the 
height z with a flux given by: 



F,{z) = -F,rr e-^^/''" , (7) 



where the optical depth Ts is: 



X*ppdz . (8) 

and the minus sign reflects the fact that this flux is going inward. The corresponding mean intensity is 
given by: 

J.(z) = ^e--A- , (9) 

The fraction of the incident radiation which is scattered is s = ap/ Xp and the fraction absorbed is 
K*p/x*p = I — s — a. For simplicity, it is assumed that a and s are constant in the region where the stellar 
radiation is mainly absorbed (both fractions are evaluated at T and P such that Tg = 2/3). This assumption 
will be checked a posteriori. 

Following CPMD, the zero-order moment of the transfer equation for the diffuse radiation field can be 
written as: 



. J Firr — -r, //In 

— — — a4:TrJs — s e , 

dTs Mo 



(10) 
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where Fs and Js are the flux and mean intensity of the diffuse stellar radiation field. 

The first moment of the transfer equation, assuming the Eddington approximation (i.e., that the diffuse 



field is isotropic), is given by (Mihalas 1978): 



From equations ( [lO| ) and ([ll|), a second order differential equation for Js can be constructed. Given 
the boundary condition corresponding to an isotropic radiation field, i.e. Js = Fs/2'k, and assuming that s 
and a are independent on Ts, the diffuse radiation field is described by: 



4^[l + (2ff/3)](l-ffVg)' 4vr(l-.gVg)' ' ^ ^ 



3[l + (25/3)](l-5V§)' +(l-5Vg)' ' 



where g = \/3a. 

Finally, from equations (0) and (|lo|), the disk heating rate per unit volume due to stellar radiation can 
be written as: 



* 

irr — XpP 



47raJ« + (l-s) — e-^^/'^« 



= 47rKj,p[J^ + J,] . (14) 



The irradiation flux Fi^r is calculated as described by KH, with the disk height and shape calculated 
self-consistently, under the assumption of gas and dust well mixed and thermally coupled. 



4. Energy transport 

In different regions of a given disk, different mechanisms could be responsible of the energy transport. 
With the aim of making a reliable model to cover a wide range of disk and stellar parameters, we have 
considered that the energy in the disk can be transported by radiation, convection and turbulence. In the 
following we describe each of these mechanisms. 



4.1. Radiative transport 



The disk's own radiation is characteristic of the disk local temperature, so it corresponds to lower 
frequencies respect to those characteristic of the stellar radiation. The radiative transport at these 
frequencies is described through the first two moments of the transfer equation integrated in frequency, 
using the Eddington approximation to close the system ( Mihalas 197S ). Instead of mean opacities averaged 
taking as weighting functions the frequency dependent mean intensity and flux (kj and xf, in Mihalas 
notation), the Planck and Rosseland mean opacities {np and xr) have been used. These mean opacities are 
computed with a consistent set of monochromatic opacities (see CPMD and D'Alessio 1996| , for details) . 
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Thc first two moments of the transfer equation in a frequency range characteristic of the disk physical 
conditions, can be written as: 



—— = Ad - 1 d = Annp p 
dz 



Jd 



(15) 



dJd ^ Fd 

dz ' An 



(16) 



In equation (|15|), represents the radiative cooHng rate per unit volume, and Ta, the heating rate per 
unit volume due to disk radiation (i.e. in the disk frequency range). 

Finally, the energy per unit area per time transported by radiation is given by the direct stellar flux, 
the diffuse stellar flux and the disk radiative flux, that is. 



(17) 



4.2. Turbulent flux 

The turbulent elements responsible for the disk viscosity are also transporting energy in the vertical 



direction (Riidiger, Elstner, fc Tschape 1988). Here we have assumed that the turbulence in the disk is 



not generated by convection but by other mechanism (e.g. instabilities associated with magnetic fields, see 



Balbus fc Hawley 199l| , 1992; [Hawley fc Balbus 199l| , 1992). So convection and turbulent fluxes are treated 



as different energy transport mechanisms, which can be both present in a given region of the disk. 

An adiabatic turbulent element moving from a hot region to a colder one in a superadiabatic medium 
has, through its trajectory, an excess of energy with respect to the surrounding medium. It loses part of its 
energy doing work on the surrounding gas and the rest of the energy is liberated when it finally dissolves, 
mixing with the medium. If the medium is subadiabatic, the rising adiabatic elements have less energy 
with respect to the ambient gas and cools the upper and cooler layers. Assuming there is no net vertical 
mean motion in the disk, there are the same number of elements rising and falling at each height, so the 
net energy interchanged between the turbulent elements and the medium along their path is zero. Only 
the excess or deficit of energy of the elements and the medium, just before dissolving, contributes to the 
net flux. If the turbulent elements are adiabatic, the turbulent energy flux is proportional to the entropy 
gradient, and following Riidiger et al. (1988), can be written as: 



Fturb — —pT 



Vt dS _ a {P + Prad) 



Pr dz 



nPr 



V 



(18) 



where S is the specific entropy, V is the gradient [d InT/d InP) of the medium, V^i is the adiabatic gradient 
(i.e. V evaluated at constant entropy) and Prad is the radiation pressure. The Prandtl number Pr is given 
by the ratio of the efficiency of momentum and energy transport by turbulent elements, i.e., Pr ~ i^t/xtj 
where vt is the turbulent viscosity coefficient and Xt is the turbulent conductivity coefficient. Finally, gz is 
the z-component of the stellar gravity, given by: 



- 8 - 



[1 + (z/i?)2]3/2 ■ (19) 



4.3. Convection 



The energy transport by convection is included through the mixing length theory (e.g. Cox & Giuli 



196S ) , with a convective efficiency computed assuming that the eddies loose energy by turbulent flux and 



by radiation. The radiative losses are calculated estimating the optical depth of the elements (Mihalas 



197S). The temperature gradient of the medium is calculated given the total flux of energy that has to be 
transported at each point. In order to calculate the energy transport by convection in both optically thin 
and thick convective regions of the disk and also to take into account that part of the total energy flux F is 
transported by turbulent conduction, instead of using the usual definition of radiative gradient based 
on the diffusion approximation, we introduce the gradient \7rc- This is the gradient (d InT/d InP) of the 
medium when convection is absent and can be written as: 

S/rc^^-j^ + Va , (20) 

Aturb 

where Aturb = ctP dz/^Pr^ A and we have used equation ([T^). 

It is assumed the disk is unstable to convection in those regions where Vrc > ^A ■ When this condition 
holds the true gradient of the medium is given by: 

V = (l-C) Vac + CVa , (21) 
where C is a convective efficiency and depends on the mixing length A, a free parameter of this convection 



theory (see D'Alessio 1996 for details). On the other hand, when V^c < '^a the medium is stable against 



convection and energy is transported only by radiation and turbulent flux. In these regions V — Vi^c■ 

The quantity V contains the information of the efficiency of the different energy transport mechanisms 
at each point in the disk. Given V, the differential equation for the kinetic temperature can be written as: 

f - -V^^^.P , (22) 

which describes the temperature structure required by the energy transport mechanisms and the assumption 
of steady state. 

The adiabatic gradient and other thermodynamical quantities used in this work have been 
calculated in the way described by Vardya (1965) and Mihalas (1967). 



5. Hydrostatic equilibrium 

The disk is assumed in vertical hydrostatic equilibrium, 

dP dPrad .„„x 

-1- =-P9z ^ — , (23) 

az dz 
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where the radiation pressure Prad is included because it could be important for disks with very high 
accretion rates or in low density regions, like the upper atmosphere. The radiation pressure gradient is 
given by (see Mihalas 197^ , pag. 170): 



dPrad 

dz 



xMv , (24) 



where Xv is the monochromatic opacity coefficient, is the monochromatic radiative flux and c is the 
speed of light. Using the mean opacities and the radiative energy fluxes defined in sections 3.3 and 4.1, the 
radiation pressure gradient can be approximated as: 

^ « ^ [XRFd + X*r{Fs + F,)] = ^ [xRiFrad " " F.) + X*r{Fs ^Fi)] , (25) 
az c c 

where XR is the Rosseland mean opacity, Xr is a Rosseland type mean opacity calculated using the 
derivative of the Planck function evaluated at the stellar temperature as the weighting function. 



6. Equations 

From the previous discussion, the set of differential equations that describe the disk vertical structure 

is: 

7— = Ad — Frf — Tirr , (26) 

az 

dJ 3 

--^^~^XR{P.T)p{Frad~F,-F,) , (27) 
dF 

'j — ^vis ~1~ Tion 5 (^^) 

az 

^ = -~V{Frad. F, T,P)^g,p , (29) 

^ =-p9.-^ ■ (30) 
az dz 



These equations are solved subject to the following boundary conditions. 

Pg=Poo, (31) 

where Pqo is a fixed small and arbitrary value of the gas pressure at 2 = z^ (we adopt P^o = 10^^ din cm^^). 
The total energy flux is: 
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The net flux produced by viscous dissipation at z — z^^ is F^is{zoa)i given by equation and the net 
flux produced by energetic particles ionization at z = z^o is Fion{zoo) from equation (||). The turbulent and 
convective fluxes are zero, since by definition the eddies cannot get out from the disk surface. The stellar 
flux going into the disk at Zoo is Fi^r- This flux is reprocessed by the disk and emerges from its surface at 
a frequency range characteristic of its photospheric temperature. From the steady state assumption, the 
incident stellar flux is equal to the emergent stellar reprocessed flux, both integrated in frequency. Thus, the 
net flux associated with the stellar irradiation is zero at Zqo and it does not contribute to equation (p^). All 
the energy flux emerging from the surface of the disk is transported by radiation, i.e. Frad{zoo) — F{zoo)- 
Finally, 



JdiZoc) — Joe — o \Fvis ^ F-ion Firr ^sMz^: (^'^) 

where the irradiation flux Firr is the stellar flux intercepted by the disk surface. The mean intensity Jq^ 
at the disk frequency evaluated at the disk surface, is calculated with the two stream approximation, 
subtracting the contribution of the diffuse radiation field to the total stellar flux because it is scattered in 
the stellar frequency range. 

At z = : 

F ^Frad^ 0, (34) 

i.e., at the disk mid-plane {z — 0) all the energy fluxes are zero because of the reflection symmetry. 

This is a two boundary problem. The disk heigh Zao is an unknown boundary and has to be 
determined. Using a 4*'' order Runge-Kutta, we solve a reduced set of equations corresponding to the 



diffusion approximation. This result is used as an initial guess for a relaxation method (Press et al. 1989) 
to integrate the full set of equations discussed above, taking z^o as an eigenvalue. 

We assume in this work that the disk is not surrounded by any substantial external material (i.e. an 
envelope or a wind). With this assumption, we choose a small value for Poo such that any material above 
it has no effect on the vertical structure; with this value of Poo, we find that the outer disk regions flare. 
However, we also find that for higher values of Poo, there can be an outer region where the disk height has 
a maximum and then decreases with P, so that the slope of the disk surface becomes negative and the the 
computed irradiation fiux becomes zero. But to have higher values of Poo, the disk should be in pressure 
equilibrium with some surrounding material. In this case, this material would be heated by the central star 
and could transport scattered and reprocessed stellar radiation to the deeper disk regions, which cannot 
"see" the star directly (Natta 1993; D'Alessio, Calvet, & Hartmann 1997). If this quasi-isotropic heating 
was included, the outer disk regions would flare (D'Alessio 1996; D'Alessio, Calvet, & Hartmann 1997). 
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Results and discussion 



For a central star of given Af» , i?, and , a disk mass accretion rate M and a viscosity parameter 
a, the density and temperature structure can be calculated by integrating equations ( p6| ) to (30). In this 
paper, we have chosen, as an example, a set of parameters that are typical of T Tauri stars in Taurus: 
M = 10-^ Mpy r"^ (|Valenti, Basri fc Johns 199^ ; [Gullbring et al. 1997| ), = 4000 K, M, = 0.5 Mq, 
and — 2 ( Kenyon fc Hartmann 1995 ). The assumed mass accretion rate is lower than the critical 
mass accretion rate that leads to thermal instability (see Kawazoe fc Mineshigc 199^ , Bell & Lin 1994), 
justifying the steady state assumption of the present model. The last parameter is the viscosity coefficient 
a. The value of a is crucial for modelling of the vertical structure of accretion disks, but is not known 
observationally. For the model discussed in this paper we chose a = 0.01. This value is consistent with 
angular momentum transport by turbulence initiated by the BH magnetohydrodynamical instability 



(Balbus fc Hawley 199l|, |Hawley fc Balbus 199l|jBalbus fc Hawley 1992| , [Hawley fc Balbus 1992i |Gammk 



1996| ). Also, as shown by D'Alessio (1996), T Tauri's observed SEDs at mm wavelengths can be fitted by 
irradiated disk models with different combinations of A/ and a. In particular, for A/ = 10^^ MQyr~^, 
the viscosity parameter which leads to a flux and slope around A = 2.7 mm consistent with observations 
(Dutrey et al. 1996) is a = 0.01. In two forthcoming papers, we will present (1) a grid of models to 
illustrate the effect of different parameters on the SEDs, and (2) a detailed comparison between models and 
observations, to determine the range of parameter space that yields the best fit to observed SEDs. 



7.1. Optical Depths 

Figure||shows the radial distribution of the Rosseland Tij, the Planck rp and the "stellar" Ts/^o{R, Zoo) 
mean total optical depths, integrated from the surface to the midplane of the disk. In this model rp > 1 
for i? < 20 AU. For larger radii the disk becomes optically thin, decreasing Tp with radius. The disk is 
optically thick to the stellar radiation (Ts/fJ.o » 1) although it can be optically thin to its own radiation 
(rp < 1). This is a consequence of the shorter wavelength characterizing the stellar radiation (CPMD) 
and the fact that it penetrates in a slanted angle (CPMD, Malbet & Bertout 1991). The Rosseland optical 
depth decreases for small radii (i? < 0.1 AU) because the sublimation of dust, which is the main opacity 
source in the disk. At these radii the disk temperature becomes similar to the sublimation temperature, 
given approximately by Tsub ~ 1800 — 2000 K. 

The difference between the Planck and Rosseland mean optical depths reflects the strong frequency 
dependence of the monochromatic opacities used to calculate the means. While rp controls the radiative 
transport in optically thick regions, rp controls the radiative cooling in optically thin regions. Because 
radiation is the most important mechanism in transporting energy inside the disk (see below) the 
characteristic temperatures are highly dependent on the radial distribution of rp and rp 



7.2. Characteristic Temperatures 

Figure || shows the temperature at the disk midplane Tc and the photospheric temperature, (where 
rp = 2/3) Tphot- For small radii {R <^ 0.5 AU in this model) flaring is not important and Tphot R~^^^ as 
in an irradiated flat disk. For larger radii, the disk curvature becomes important, increasing the amount of 
stellar flux intercepted and flattening the distribution of photospheric temperature, giving Tphot R^^^^. 
This result is consistent with the temperature distribution found by Miyake fc Nakawa (1995), for an 
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optically thick flared disk, with gas and dust well mixed. We only show the photosphcric temperature for 
regions where tr < 2/3; outside, it is not defined. A simple analytical approximation gives Tphot R~^^'^ 
for the region where the disk flaring becomes important, assuming that the stellar radiation is deposited 



at a height Zg c>c where H is the gas pressure scale height (see Canto, D'Alessio fc Lizano 1995 ). In 
this paper, the disk surface is assumed at constant pressure instead of proportional to the scale height, 
which results in a different radial dependence of the irradiation flux. We will show below that Zqo is not 
proportional to H. 

For R <j2 AU, viscous dissipation is the main energy source of the regions near the disk midplane. At 
larger radii, Tc becomes higher than the central temperature of a disk heated only by viscous dissipation, 
reflecting the importance of irradiation as a heating source. The contrast between Tc and Tphot depends 
on Tn, since this optical depth controls the vertical temperature gradient required to transport the energy 
released at every height. 

Figure || also shows that the central temperature distribution becomes flat around Tc ^ 1800 — 2000 K 
{R <^ 0.1 AU) where the Rosseland mean optical depth has a minimum (see Figure |^). The dust is the 
dominant opacity source at temperatures lower than the sublimation temperature, which depends on the 
type of grain and the local density. In the regions where dust is destroyed the Rosseland mean opacity XR 
decreases. At a given annulus, the optical depth increases towards the midplane (i.e., with a decreasing z), 
until a depth is reached where the dust is destroyed. From this height to the midplane the optical depth 
remains constant. In the region where dust is destroyed, the temperature is almost constant with height, 



being Tc ~ Tsub (see D'Alessio 1996). If the density were high enough to increase the optical depth towards 
the midplane (if there is a heating source at the disk midplane, e.g., viscous dissipation) then Tc > Tsub- 
This is the case of disks with a lower a or a higher M than the model presented in this paper. 

As can be seen in Figure ^ the Rosseland optical depth decreases with radius, so the gradient 
of temperature required to transport the energy produced in the disk also decreases with radius. The 
central temperature, for the optically thick annuli where the diffusion approximation holds, scales as 
Tc ~ Tphot (3t/j/4)^/^. For i? ^ 10 AU, where r/j <^ 10, the photosphcric and the central temperatures 
becomes very similar, i.e., the disk interior is almost vertically isothermal. In this outer region, the 
temperature is given approximately by T w O.SA{Firr / (Jr)^^'^ , as follows from equation (^), with 
J = Joo = const, T = const, assuming rp >> 1, Firr >> Fms, Fion, and neglecting the diffuse radiation 
fleld. 

At a given annulus, the density decreases with height and thus the optically thick annuli have upper 
optically thin layers, i.e., a disk atmosphere. Figure^ shows the vertical distribution of temperature of the 
reference model at different radii. For all the disk there is a temperature inversion produced by the stellar 
irradiation (see CPMD). This temperature inversion is due to two effects, previously explored by CPMD: (i) 
the stellar radiation penetrates the disk along a slanted direction so the optical depth along its trajectory is 



larger than the case of a penetration perpendicular to the disk surface (see also Malbet & Bertout 1991), 
and (ii) the opacity of the disk material (dominated by dust) to the stellar radiation, is larger than the 
disk opacity to its own radiation (due to the shorter frequencies of the former). Malbet & Bertout (1991) 
considered only the flrst effect, but the second effect is dominant as can be seen in Figure |[ 

The main difference between CPMD's treatment and the model presented in this paper is that we 
calculate the disk optical depth self-consistently with the temperature and density vertical structure, while 
CPMD assume the disk interior is optically thick and that there is no viscous dissipation at the disk 
atmosphere, where energy is only transported. With the same set of equations, we calculate the optically 
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thin atmosphere and the optically thick interior structures, without any assumption a priori about the 
optical depth of the different regions. Annuli situated at large distances from the star, can be completely 
optically thin and the assumptions made by CPMD are not valid, while with the method described in 
this paper, the structure of the optically thin outer disk can be calculated. Thus, we have found that 
these regions have a temperature inversion, being the surface temperature higher than the photospheric 
temperature. The annuli with low Rosseland mean optical depth, are almost vertically isothermal between 
the midplane and the photosphere, but has a large contrast in temperature between the photosphere and 
the disk surface. 

Figure ^ also shows that the approximation found by CPMD, with a correction introduced by the use of 
the Planck mean opacity in equation (p^), is a good description of the temperature in the disk atmosphere. 
The disk upper atmosphere, with a temperature higher than the photospheric temperature, is like the 
"super-heated" layer proposed by Chiang & Goldreich (1997), where the stellar radiation penetrating the 
disk is deposited. The main difference with their treatment and this paper is that here the disk vertical 
structure is calculated in detail. 

Figure |5| shows the surface temperature of the disk Tq, i.e., the temperature at Zoc- From R 0.5 AU 
to 80 AU, To is given approximately by Tq ^ R~^-^^, and for larger radii we found Tq ~ R~^^^, which is 
consistent with the expected temperature of optically thin dust (with an absorption coefhcient kx cx: 1/A) 



heated by stellar radiation (Spitzer 1968 , pag. 194; see also Chiang fc Goldreich 1997) 



7.3. Characteristic Heights 

There is not a unique definition of disk height. Figure |^ shows the disk photospheric height Zphot, the 
scale height of the gas at the midplane temperature H, the height at which the stellar radiation is deposited 
Zs, and the height at which P = Poo given by Zoo- There are clear differences between these heights which 
can be understood as follows. 

The pressure scale height of the gas at the central temperature is given by: 

H ^ c,(T,) 
R Rn{R) 

If Tc/ fj, decreases with R slower than 1/R, then the gas pressure scale height increases with R. The 
higher the central temperature, the larger the pressure scale height. Therefore stellar irradiation increases 
this scale height with respect to a purely viscous disk. For R ^ 6 AU, H cx i?^/* as it is expected from 
equation ( p5| ) when Tc oc i?^''-^. 

The photospheric height Zphot is defined as the height where r/? = 2/3, if the disk is optically thick. For 
i? J> 30 AU the disk becomes optically thin and there is no photosphere. The ratio Zphot/H w 2.5 between 
0.1 and 20 AU. 

The absorption height z^, given by the height at which Ts/ho = 1, is the height where the largest 
fraction of stellar radiation is absorbed. Figure |6| shows that Zs is larger than Zphot at every radius of the 
disk, reflecting the fact that the disk becomes optically thick to the stellar radiation above the depth at 
which it becomes optically thick to its own radiation. For i? > 20 AU the disk is optically thick to the 
stellar radiation and optically thin to its own radiation. At large distances from the star, z^ is close to the 



kTcR 
GM^fi(Tc, Pc)mH 



(35) 
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maximum height of the disk model Zao- 

From the vertical structure of this reference model, we find that z^/ H decreases from 4.5 to 2.5, when 
R changes from 10 to 300 AU. The ratio Zg/H is assumed constant in models where there is no vertical 



structure calculation (e.g. KH, Canto et al. 1995, Chiang & Goldreich 1997) leading to Tphot oc R ^''^ 
instead of Tphot oc 



The physical model considered in this paper is different from that described by p3ell et al. 1997 , where 



the disk is illuminated by a uniform radiation field with different temperatures (Tji — 10, 20andl00 K). 
In that case, the relevant height is the vertical optical depth at the characteristic wavelength range of the 
incident radiation. Bell et al. find for the lowest radiation temperatures, the height decreases at large radii 
so that the outer disk regions cannot see the star. In the model studied in this paper, he radiation is coming 
from the central star (T* = 4000 K), with a large incident angle respect to the normal to the disk surface 
and also at shorter wavelengths. Both properties combine to produce a flared disk, with an absorption 
height increasing with radius. Therefore, when stellar irradiation is included in a self-consistent way, taking 
into account its characteristic wavelength and direction, the outer disk can be heated by the central star. 

We can see that all the characteristic heights plotted in Figure |^ are less than R, supporting the thin 
disk approximation, used to decouple the vertical and radial structures. 



7.4. Surface density 

Given the disk density as a function of R and z, the disk surface density radial distribution Y.oo{R) is 
calculated using equation (^). Figure |^ shows Soo(7?) for an irradiated and a non irradiated disk with the 
same M , a and central star. From the equation of conservation of angular momentum in a steady a— disk, 
the surface density can be written as: 



M 



Sir < Vt > 



R 



1/2- 



(36) 



where < ft >= ^'^^vtiz, R) p{z, R)dz/T,ao is a mean viscosity coefficient. Approximating < > by the 
viscosity coefficient evaluated at the disk midplane, the surface density can be written as: 



M Vt^niH f R* 
Sira kTc \ R 



1/2. 



(37) 



Figure |^ shows this expression for Eqo compared with the surface density calculated numerically from 
p(i?, z). The approximation is a good description of the numerical result, which can be used as a check of 
the method used to integrate the vertical structure equations. From i? = 2 to 100 AU, Sqo ~ as it is 
expected from equation ( ^7|) when Tc ~ Then, a steady a-disk has a surface density distribution 

flatter than the nominal form Sqo ~ {R/ R^,)~^^^ , usually found in the literature (e.g. [Bcckwith et al 



1990| ; putrey et al. 1996| ; [Miyake fc Nakagwa 1995| , [Chiang fc Goldreich 1997| ). Figure g also shows this 
power law surface density for a disk with the same total mass (assuming R^ = 100 AU) than our reference 
model. The a irradiated disk has denser outer regions than the nominal model. This implies that the 
brightness of the outer disk is higher (for a fixed temperature distribution and opacity) at submm and mm 
wavelengths, where the disk becomes optically thin. This result has implications on the physical properties 
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of disks as inferred from observations at long wavelengths, which will be discussed in paper III. 



7.5. Disk mass 



The mass of a disk model is given by: 



IB,, 



(38) 



where Ri is the inner radius of the disk, which can be R■^, if the disk extends inward to the surface of the 



star, or the radius of an inner hole when a magnetosphere is present (e.g. Gosh & Lamb 197£), and 
Rd is the outer radius of the disk. Due to the increase in surface area, Md is determined by the mass of 
the outer regions, in spite of the fact that the surface density decreases with i?, as discussed above. From 
equation (^ ) we can see that unless the unlikely case of the central temperature increasing with i?, the 
area of the annuli increases with R more rapidly than the decrease in Sqo. This result implies that the 
disk mass is not sensitive to the value of Ri (if Ri << Rd) and that the heating mechanisms other than 
viscosity, affect strongly the mass of a disk model with fixed values of AI and a. 

Figure || shows the mass of the disk as a function of Rd- At Rd ~ 100 AU the disk mass is 
Aid ~ 0.017 Mq, larger than the Minimum Mass of the Solar nebula, ~ 0.01 Mq, but smaller than the 
central star mass A/* = 0.5 Mq. Therefore, it seems to be a good approximation to neglect the disk self 
gravity with respect to the stellar gravity, as was assumed when the angular velocity of the disk was taken 
as the Keplerian angular velocity and, also, in the hydrostatic equilibrium equation. 

Based on the result that the disk mass is determined by its outer regions, a rough estimation of Md 
can be constructed assuming that the disk is vertically isothermal. Taking the midplane temperature as a 
power law, Tc ~ Ti{R/ Ri)~^ for R> Ri with Ti = Tc(i?i), and equations ( ^7| ) and (|3^, the disk mass can 
be written as: 



Md « M„„er(i?i) + ^^VGM.Ri — h + l/2)-i 
6k1 1 a 



Rd 
R'i 



7+1/2 



(39) 



where Minner{Ri) is the disk mass corresponding to regions with R < Ri, which we find is 



Mi, 



^{Ri) « Md{Rd) if Rd » Ri- For the reference model we find i?i « 456 i?*, Ti « 58 K 



and 7 ~ 0.5. The mass Minner{Ri) = 5 x 10 ^ M©, and Aid oc Rd- Figure || shows the approximation 



described by equation (39) compared with the numerical result. 



7.6. Characteristic timescales 

The timelife of a disk with M = 10""^ Moyr"! and Aid = 0.017 Mq is Md/AI « 1.7 x lO*' yrs. The 
characteristic time of the matter diffusion due to viscous torques tyis, is: 

R R? R?Vl 

tvis « 1 r ~ = 7777T (40) 

\ur\ <vt> aCs{Tc) 
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where is the radial accretion speed. 

The viscous time evaluated at Rj^ is an estimation of the lifetime of a disk which receives little mass 
from the surrounding cloud. Figure ^ shows that the reference model at 100 AU has tyis « 2 x 10^ yrs, also 
consistent with the disk mean age ~ 10® yrs, inferred from the position of T Tauri stars in the HR diagram 
(Strom et al. 1995). 

From equation ( ^ ) we can see that the larger the central temperature the shorter the lifetime of the 
disk, for a given value of a. Therefore, the stellar irradiation, which is the most important outer disk 
heating mechanism, can not be neglected in a lifetime estimation of a given disk model. 

7.7. Energy transport 

Figure ^ shows the fraction of the total energy flux transported by radiation, convection and 
turbulence, as a function of height, at various distances from the central star. Also the position of the gas 
pressure scale height H is shown. 

A flux with a positive sign is going in the +z— direction, and with a negative sign, is going in the 
—z— direction. In the convective regions (see the panel corresponding to i? = 0.093 AU of Figure |l^) the 
turbulent flux is positive, thus it is transporting energy from inside to outside (form inner regions to upper 
regions) . 

Convection only transports an appreciable amount of energy near the midplane around 
i? = 10 i?* ~ 0.01 AU. In the radiative zones, the turbulent flux becomes negative, as can be 
seen from equation ([l8| ) when V < V^i. In these radiative or non-convective regions, the turbulent flux 
transports energy from the upper layers, where temperature and pressure are low, to the inner layers, which 
are hotter and have a higher pressure. When a turbulent adiabatic eddie moves in a subadiabatic medium 
from a cold to a hot region, it will have a higher temperature respect to the surrounding medium. On the 
other hand, if this adiabatic eddie moves from a hot to a cold region, it will have a lower temperature 
respect to the surrounding medium. The work made by the medium on the element or by it on the medium, 
is responsible for this behavior. Therefore, the presence of the turbulent flux increases the radiative flux, as 
can be seen in Figure |l^. At every annuli, we can see that radiation is the most important energy transport 
mechanism inside the disk, enhanced in the non-convective regions because the presence of the turbulent 
flux. Thus, radiation is controlling the disk temperature structure. 

Note that less than 70 % of the total flux is produced below the gas scale height. On the other 
hand, above the disk photosphere, the total flux is almost constant (see the panels corresponding to 
R = 0.093, 0.93, 9.3 AU from Figure |l^). So, it seems to be a good approximation to neglect the 
flux produced there by viscous dissipation and ionization due to energetic particles. Then, in the disk 
atmosphere, the energy is being only transported. 

7.8. Disk gravitational stability 

Other interesting physical property that can be calculated from a detailed disk model is the gravitational 
stability parameter introduced by Toomre (1964), which quantifies the fiuid reaction against axisymmetric 
gravitational perturbations. The Toomre parameter can be written as: 
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such that, if Qt > 1 the fluid is stable. 

Figure |l^ shows Qt as a function of R for the reference model; the disk is stable for the range of radii 
used to calculate the disk structure {R < 337 AU). It is interesting to note that if the irradiation flux is 
higher Qt increases due to the increase in disk temperature and the decrease in surface density. Thus, 
irradiation tends to stabilize the disk against gravitational perturbations. 



7.9. Ionization structure 

Given the vertical structure, it is interesting to evaluate where the ionization fraction (controlled by 



energetic particles and thermal ionization) is not enough to sustain the BH instability (Balbus & Hawley 



1991| ). Gammie (1996) has called this region "the disk dead zone", where turbulent viscous dissipation 



cannot exist if the source of turbulence is the mentioned magnetic instability. 

Following Gammie (1996), the ionization fraction which corresponds to a magnetic Reynolds number 
RcM = VaH/t] = 1, is given by: 

where Va is the Alfven speed, approximated as Va ~ a^/^Cs(T), rj is the resistivity, given by 
1] = 6.5 X lO'^a;"^ cm^ s^^, x — n^/riH is the ionization fraction, and is the electronic density. 

On the other hand, from Stepinski (1992), the local ionization fraction due to energetic particles, 
considering recombination upon grain surfaces and ion reactions, can be written as: 




= 5.2 X 10-« ( ^/ 1 + ,Q-, '.riCra,i''AlHCcose^^^') ^^^^ 

where rgr is the grain typical size (we adopt Vgr = 0.1 /im). The cosmic rays ionization rate decreases 
exponentially with the surface density, with a characteristic attenuation density given by ^ 100 gcm^^. The 
thermal ionization fraction xth is calculated assuming LTE, and using the procedure described by Mihalas 
(1967). We approximate the ionization fraction in the disk as a; ~ max{xth,Xep)- 

Figure ^ shows a contour of x/xi — 1. In the region where x < xi, the disk material is decoupled 
from the magnetic field, Re^ < 1 and the BH instability cannot be the source of turbulence. This is the 
"dead zone", where there can be no viscous dissipation, and mass transport. In spite of the disk heating 
due to stellar irradiation, the thermal ionization is not enough in this region to make Rcm < 1. If there is 
no other source of turbulence, the region between 0.2 and 4 AU cannot be treated with a constant M and 
a. Outside this region, the magnetic Reynolds number is Rcm > ^ and the disk turbulence can be initiated 
by the BH instability. 

As can be seen in Figure the region where x < xi extends above the disk scale height, but it remains 
below the photosphere and the height where the stellar radiation is deposited. For R « 2.5 AU, the surface 



- 18 - 



density in the dead zone is « 90 % of the total surface density at this radius. At i? ^ 2 AU the surface 
density drops below the attenuation scale for cosmic rays, ~ 100 g cni~^, but the ionization fraction implied 
by equation (^) is smaller than the upper limit used by Gammie (1996), and Xep remains less than xi 
between 2 and 4 AU. 



7.10. Spectral energy distribution 

For the reference model we have calculated the SED, assuming the disk is pole-on. The maximum and 
minimum radii adopted to calculate the SED are Rd = 300 AU and Rhoie = 1 -R*, respectively. For the 
millimeter-wave dust opacity, we use a power law given by: 

^^^ = ^^"(20^) 

where kq — 0.07 cm^g~^ is obtained assuming that at wavelengths shorter than 200 /im the dust opacity is 
given by Draine & Lee (1984). The exponent /? is taken as 1. This opacity is, at A = 1.3 mm a factor of 
2.14 smaller than the opacity used by Beckwith et al. (1990). 

Figure ^ shows the predicted SED compared with the SED of AA Tau, a typical Classical T Tauri 
star and the agreement is good for A > 2 jim. We want to emphasize that this is not a best fit model for AA 
Tau, is only a typical model. In paper III we will use observational restrictions as the mass accretion rate 
and central star properties inferred from optical observations to restrict the family of models for different 
particular objects. Nevertheless, the comparison between the SED of a typical model with the SED of a 
typical CTTS show that the model predictions agree well with observations. 



8. Effect of stellar irradiation on the disk vertical structure 



Disks can be irradiated by external sources, such as an infalling dense envelope (Butner, Natta & 



Evans 1994; D'Alessio, Calvet & Hartmann 1997) or a tenuous dusty envelope, remnant of the infalling 



envelope or associated with a disk or a stellar wind (Natta 1993), or the central star (Adams & Shu 1986; 



KH; CPMD and this paper). Different sources can be important in different evolutive phases of a disk, and 
in this paper we explore the case of a disk irradiated by the central star. A non-irradiated model is used to 
quantify the effects of stellar irradiation on the disk structure and physical properties, although it does not 
represent a realistic disk for the case treated here (with M ~ 10^^ M0 yr^^). 

Irradiation is the main heating source of the disk. Also the heating due to ionization by energetic 
particles becomes larger than the viscous heating for R > 200 AU, but it is negligible respect to the stellar 
irradiation. This behavior of the flux is reflected in the photospheric temperature distribution, plotted in 
Figure |^. For R ^ 0.5 AU, Tphot is similar to that of an irradiated "flat disk" . For larger distances, becomes 
flatter than the temperature distribution of a viscous disk or a flat irradiated disk, Tyis ~ Tfiat oc i?^'^/^. As 
shown in Figure |l^, the irradiated disk has a larger emergent flux and a flatter SED at every wavelength, 
than the non-irradiated disk. For R <^ 2 AU, Tc is the same for the irradiated and non-irradiated models, 
because the viscous dissipation is the dominant heating source of the regions close to the disk midplane. At 
larger distances, the stellar radiation can penetrates the disk and it becomes the most important heating 
source at every height. 
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Figure |^ shows the surface density of the irradiated and non-irradiated disk model, with the same M 
and a. SteUar irradiation decreases the disk surface density for i? ^ 1 AU. As can be seen from equation 
(^), the change in disk temperature caused by the irradiation substantially increases the viscosity, and so, 
for a fixed mass accretion rate, the surface density drops correspondingly. The mass of the disk, determined 
by its surface density at large radii, also decreases with irradiation. As can be seen in Figure ||, the mass of 
a pure viscous disk with M = 10"^ Mq yr-\ a = 0.01 and Rd = 100 AU, is Md = 0.063 Mq, approximately 



a factor of 4 larger than an irradiated disk with the same M and a. From equation (|39D, a non irradiated 
disk has Md ~ PlJ'^ ^ assuming 7 w 3/4 in the outer regions, so its mass increases with the disk radius more 
rapidly than the mass of an irradiated disk Md oc Rd- 

The increase in temperature produced by irradiation decreases the viscous time of the disk, as can 
be seen in Figure 0. The non-irradiated disk has « 10^ yrs, an order of magnitude larger than the 



estimated mean timelife of disks (Strom, Edwards, & Skrutski 1993) 



Another property for which the effect of irradiation is important is the Toomre instability parameter. 
Figure |ll| shows that the non-irradiated disk becomes unstable for R ^ 28 AU. The higher temperature 
and lower surface density of an irradiated disk, with respect to a pure viscous disk with the same a and M, 
stabilizes the disk against gravitational perturbations for (at least) R <; 337 AU. 



9. Effect of irradiation from an accretion ring on the disk structure 



The magnetospheric model has become widely accepted (e.g. Konigl 1991; Calvet & Hartmann 1992; 



3hu et al. 1994 ). It is expected that the presence of a magnetosphere modifies the disk model discussed 
in this paper. Two important effects on the disk structure and emission are: (1) to truncate the disk, 
such that the resulting inner hole decreases its near IR emission and (2) to heat the disk, changing its 
temperature and surface density distributions. The disk heating is produced by the the interaction with 
the stellar magnetic field ( Kenyon, Yi fc Hartmann 1996| ) and also, by radiation from an accretion shock 



at high stellar latitudes. For low mass stars, this radiation is characterized by photons of much shorter 
wavelength than those of the stellar radiation or the local disk emission. Kenyon et al. (1996) conclude 
that the effect of this ring radiation is small, but they do not include the increased absorption efficiency 
due to the wavelength-dependence of the opacity. In this section, we study the eff'ect of the accretion ring 
radiation on the vertical structure of the reference model. The dependence of this additional heating source 
on disk and stellar parameters and its effect on physical (mass, stability, etc.) and observable (near IR 
colors, long wavelength fluxes, etc) properties of disks will be calculated in papers II and III. 

The calculation of the ring flux Fa intercepted by the disk and its mean incidence angle cos^^ is 
described in the Appendix. The ring is located at an angle ^0, measured from the z axis (see Figure 
and has an angular width A^o- We assume the ring has a temperature Ta = 10000 K, inferred from the UV 



excess of T Tauri stars (e.g. Bertout, Basri fc Bouvier 1988). The ring angular width can be estimated as: 



A...0^057 (,i„W- il)(^) (5I-) ru, ,45) 



where La is the ring luminosity (for a ring seen pole-on by the observer). The ring luminosity can be 
related to the mass accretion rate (see Kenyon et al. 1996), writing: 
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The disk adopted as a reference model, has AOq = 0.062°/(sin6'o)~^ {2 — R^, / Rhoie) ■ We have taken two 
values for 6*0 = 30°, 60°, and an extreme value (2 — R^,/Rhoie) ~ 2, to explore the effect of the irradiation of a 
ring on the disk structure. The corresponding ring widths arc: A^o = 0.25°, for do = 30°, and Ado = 0.14°, 
for ^0 = 60°. It is clear from equation (^) that the higher the M, the larger the ring luminosity, and the 
larger effect on the disk, at least on its upper layers or outer regions, where viscous dissipation can be less 
important. 

We preserve the same set of equations and boundary conditions given in §6, but now the irradiation 
flux Firr is substituted by the stellar irradiation flux plus the ring flux, i.e.. 



Frr ^ {I - nA)F.*.r + Fa, (47) 

where flA is the fraction of the solid angle subtended by the star, occupied by the ring (see the Appendix), 
and F*.j. is the irradiation flux from the star (see KH). 

The cosine of the angle of incidence is estimated as: 

^T4(l - ^a) + ^TXQa 
MO - <! ^ } , (48) 

mi-nA) + TXnA 



where cos^^{^a) is the averaged angle between the direction of the ring flux and the normal to the disk 
surface (see the Appendix), and cos~^(/i*) is the mean angle between the stellar radiative flux and the 
normal to the disk surface. The mean Planck mean opacity, which describe the interaction between the disk 
material and the incident radiation fleld, is written as: 



K*pT^{l - ^a) + K^T^nA 
= \ ^ \ , (49) 

n{i~nA) + TXnA 



and the other mean opacities ^md &re calculated with the same type of average. The mean "ring" 
opacities Kp, Xp and Xri ars the mean opacities calculated using the Planck function evaluated at Ta as 
the weighting function (see eq. [^). 

At every radius, the irradiation flux from the star is larger than the irradiation flux from the ring. For 
R > 0.1 AU, the stellar flux is a factor of 10 larger than the ring flux. The ring irradiation increases the 
disk photospheric temperature less than 1 %, in agreement with Kenyon et al. (1996). The temperature at 
the disk surface increases less than 6 %. 
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10. Summary and conclusions 

In this paper we discuss a model of the vertical structure of steady and thin accretion disks, irradiated 
by their central stars, that includes transport of energy by radiation, convection and turbulent flux. The 
radiative transport is treated without the diffusion approximation: The optically thin and thick regions 
of the disk are calculated with the same equations, without an a priori separation. We discuss the 
physical properties of a disk model with M — 10^^ Mq yr~^ and a — 0.01, around a low mass star with 
Af* = 0.5 Mq, i?, = 2 R© and T* = 4000 K, which represents a typical T Tauri disk. The effect of other 
input parameters on the structure and physical properties of disks will be reported in paper II. 

From this study, we conclude that the stellar irradiation is the main heating source of the outer regions 
(i? > 1 — 2 AU) in a typical T Tauri accretion disk, and has to be taken into account to calculate its 
physical properties, such as: temperature, surface density, mass, height, timelife and gravitational stability. 
For instance, the photospheric and midplane temperature distributions of an irradiated disk are flatter 
than those of a non irradiated viscous disk. The former case has a smaller surface density and, therefore, 
a smaller total mass than a non irradiated viscous disk. Moreover, irradiation tends to stabilize the disk 
against gravitational perturbations, increasing the Toomre's parameter above one. Therefore, if some of 
the disk properties are inferred from observations, they can be related to the disk parameters M and a, 
through a detailed model, that must include the stellar irradiation. 

Stellar irradiation produces a temperature inversion at the disk atmosphere, which, as shown by 
CPMD, has observational consequences. The sensitivity of the disk upper atmosphere temperature to the 
mean opacity, suggests that it is important to make a detailed "model atmosphere" of the disk, solving 
the radiative transfer equation at different wavelengths and directions. Moreover, we have assumed that 



gas an dust are well mixed and thermally coupled. HST observations of HH30 (Burrows et al. 1996) 
indicate that the dust scatters light from the central object and producing the observed image at scale 
heights H/R < 0.1, supporting the assumption that gas and dust are well mixed. Nevertheless, in order 
to evaluate both assumptions it seems necessary to calculate the disk structure considering gas and dust 
as separate ingredients with their particular dynamical and thermal behavior. This is a problem outside 
the scope of this paper, but we want to emphasize the importance of both assumptions. On one hand, 
the gas is heated by viscous dissipation while dust is heated by stellar irradiation. If the coUisional and 
radiative interaction between gas and dust is efficient enough, both can be characterized by the same 
temperature. On the other hand, if gas and dust are well mixed by turbulent motions, they have the 
same spatial distribution. In any case, given both assumptions, the height of the disk can be calculated 
self-consistently including the irradiation. A simple estimate of the difference in temperature between gas 
and dust, including only radiative and collisional heating of the gas particles by dust grains, shows that the 
assumption of a unique temperature is good below the disk photosphere. Nevertheless, the region where 
the stellar radiation is mainly absorbed has such low pressure that coupling is not very efficient. This 
can introduce an inconsistency in the atmospheric calculation discussed in the present paper, which has 
to be evaluated with a separate calculation of gas and dust temperatures. Also, as discussed by Chiang & 
Goldreich (1997), to calculate molecular line emission it is important to take into account the possibility 
that gas and dust have different temperatures in the upper atmosphere. In this sense, the present model 
corresponds to a first approximation to the atmospheric structure, which has to be evaluated with a more 
detailed treatment of the gas and dust thermal and dynamical behavior. 



The accretion ring on the stellar surface can irradiated the disk. For the reference model this ring has 
a luminosity around 0.02 — 0.04 Lq and a temperature around Ta = 10000 K (from the UV excess of T 
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Tauri stars). In this case, we found that the effect of this additional heating source on the disk structure 
and emission, is neghgible. 

Our formahsm can include any physical model that provides the viscosity v as a, function of local 
physical conditions in the disk. For the particular model discussed in this paper, we have used the a 
prescription, assuming a constant value of a throughout the disk, to evaluate v. Given the vertical structure 
and following Gammie (1996), we have found the region where the ionization fraction (including thermal 
and energetic particles ionization) corresponds to a magnetic Reynolds number of Rem < 1, for a disk 
model with constant a. In the region where Rcm < 1 the magnetic instability proposed by Balbus & 
Hawley (1991) cannot work. If the source of turbulence is this magnetic instability, then there cannot be 
viscous dissipation where Rcm < 1, i-e., for 0.2 < i? < 4 AU, and this region will represent the dead zone 
of Gammie's model. In particular, Gammie proposes a "layered accretion" throughout this dead zone to 
alleviate the problem of the non-steady behavior of accretion disks. 



Finally, the predicted SED for the reference model is in agreement with the SED of AA Tau, for 
which preliminary results from Gullbring et al. (1997) are: M = 3.8 x 10~^ Mgyr^^ = 0.53 Mq, 
i?, = 1.74 Rq and = 4048 K. Both SEDs are shown in Figure |l|. The wavelengths at which the SED is 
plotted, are thought to sample the continuum spectrum, and not to resolve lines and bands. Nevertheless, 
there are some apparent features in the irradiated disk SED, like the CO fundamental tone in emission 
around 4.3 /im and an IR band of water vapor around 100 /im. It is important to emphasize that these 
bands in emission are due to the atmospheric temperature inversion, and the assumption that gas and dust 
are thermally coupled. If this last assumption is not valid at the upper disk atmosphere, the details of the 



gas spectral features will change respect to what is shown in Figure O (see Chiang & Goldrcich 1997) 



Details of the disk SEDs, calculated for different inclination angles, will be discussed in a forthcoming paper. 



We are grateful to Javier Ballesteros, Lee Hartmann, Alejandro Raga and Salvador Curiel for helpful 
discussions. This work was supported in part by Instituto de Astronomi'a, UNAM, Mexico, DGAPA-UNAM, 
ConaCyT, NAGW 2306 and NAG5-4282. 



- 23 - 



A. Irradiation of the disk by a hot ring on the stellar surface 

In this appendix we describe the calculation of the flux intercepted by the disk, emergent from a hot 
annulus on the stellar surface. This annulus would correspond to the accretion shock between the material 
flowing through the stellar magnetosphere and the surface of the star. Figure |lj shows the geometry of this 
problem. The coordinate system and angles are defined in this figure. 

A surface element of the hot annulus is given by: 

dA = rI dfi (sin 9 cos + sin 9 sin cpy + cos 9z) (Al) 

where d^ = sin 9d9d(j> is the solid angle of the surface element as seen from the center of the star. The same 
surface element, seen from a point P on the disk surface, has an apparent area dAp, given by the projection 
of dA on d, where d = dp — f, dp is the vector describing the position of P, and r is the position vector of 
the annulus element, both given by: 

f = R^, (sin 9 cos (fix + sin 9 sin (j)y + cos 9z) (A2) 



dp^Ry + ZocZ (A3) 
Thus, the unit vector d can be written as: 

J dp — r — R^, ain 9 cos 4>x + {R— R* sin 9 sin (l))y + (zoc ~ R* cos 9) z (Ai) 
~ |4 - 7=1 ~ + i?2 + _ 2RR, sin 6* sin - 2zooR* cos 9] 1/2 ^ ' 

The annulus element, as seen from point P, has a solid angle given by: 

dAp dV,Rl [R sin sin — i?, + Zqo cos 9] 

^ ~ "iP' ~ [Rl + R^ + zl^- 2RR^ sin 6I sin - 2zooR* cos 9]^/^ ^ ' 

The mean intensity of the annulus radiation field at point P is given by: 

The flux intercepted by the disk is: 

Fa= [ [ I^^ld^pdv = / iid^p (A7) 



where Ta is the ring effective temperature and is the cosine of the angle between the vector d, which 
characterizes the direction of incidence of the radiation, and the vector n, normal to the disk surface. The 
disk normal vector can be written as: 
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. _ -{dz^/dR)y + z 



Thus, 



t . [(i? — i?, sin0siii(7!))(izoo/rf^ — Zoo + -R* / , „n 

u = — d • n = — ( A9) 

^ {l + {dzoo/dR)^y/'^[Rl+R'^ + zl,~2RR^sin9sm(l)-2zooR*cose]^/^ ^ ' 



Finally, substituting equation (A9) in (A7), the flux of the annulus intercepted by the disk can be 
calculated. The averaged value of the cosine of the angle between the incident radiation and the disk normal 
can be calculated as: 

f d ■ ndn„ 

^^ = -Jdnf 

For simplicity, we assume that the annulus is at ^ = 9q and it is thin enough to approximate the 
integrals over 9 by the integrand evaluated at 9q times 59q. For R < (i?, — z^o cos^o)/ sin^O: the annulus is 
no seen by the disk. The integral over cj) is symmetric respect to (j) = 7r/2, then we make it from 4'min to 
4>max, and multiply the result by 2. The limits of this integration are given by: 



1 R^ Zqq cos 

= sm 

R sin 9o 

if \R, 

^oqCOS^oI — |-Rsin^o|i otherwise (j)mi7i — (ind (ptnax 

= 7t/2 



>™„ = sm „ . „ (All) 

R sm t^o 
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Fig. 1. — Geometry of the disk. This plot shows schematically what is called annulus, surface, atmosphere, 
interior and midplane through the text. At a given distance R from the central star, the disk vertical 

structure refers to the disk physical properties in the z— direction. The atmosphere refers to the optically 
thin part of the disk and the interior corresponds to the optically thick part; we indicate schematically that 
the outer regions become optically thin. 
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Fig. 2. — Radial distribution of different mean optical depths, integrated from the disk midplane to its 

surface. The disk parameters arc M = 10^^ MQyr"-'^ and a = 0.01, the central star parameters arc 
M* = 0.5 Mq, i?* = 2 Rq and = 4000 K. The vertical structure was calculates from: R = 0.0127 AU 
to 373 AU. The plotted optical depths are: Rosseland mean tr (solid line), Planck mean rp (dashed line) 
and the mean optical depth at the stellar frequency range Ts / ijlo{R, Zoo) (dotted line). The horizontal line 
corresponds to r = 1. 
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Fig. 3. — Radial distribution of the midplane and photospheric Tphot temperatures for an irradiated and 
a non- irradiated disk model. The plotted temperatures are: Tc for the irradiated disk (solid line), for the 
non-irradiated disk (dotted line), Tphot for the irradiated disk (stars) and Tphot for the non- irradiated disk 
(triangles). The disk and stellar parameters are the same than Figure ||. 
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Fig. 4. — Vertical temperature distribution of an irradiated disk vs logarithm of the Rosscland mean optical 
depth (integrated from the disk surface to z). The curves represent: the temperature calculated in this 
paper (solid line), the temperature of a non- irradiated viscous disk (dashed line), the temperature given by 
CPMD (dot-dashed line), temperature calculated with the method of Malbet and Bertout (1991) (dotted 
line). Irradiation produces a temperature inversion towards the surface at all radii, the photosphere is at 
T/f = 2/3. The disk and stellar parameters are the same than Figure ^. 
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Fig. 5. — Temperature at the disk surface Tq, for an irradiated disk (solid line) and a non irradiated disk 
(dotted line). The photospheric temperature Tphot of the irradiated disk is shown as a reference (stars). The 
disk and stellar parameters are the same than Figure 0. 
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Fig. 6. — Radial distribution of different characteristic heights of an irradiated accretion disk. The plotted 
heights are: the gas scale height H evaluated at the central temperature (solid line), the disk surface height 
Zoo (dashed line), the height where a large fraction of stellar radiation is absorbed (dotted line), and the 
photospheric height Zphot (dot-dashed line), all in units of the disk radius. The disk and stellar parameters 
are the same than Figure ||. 
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Fig. 7. — Radial distribution of surface density of an irradiated (solid line) and a non irradiated (dashed 
line) disk model. For the irradiated disk, the approximation for Y^oo given by equation ( ^7| ) is also plotted 
(stars). The nominal surface density E cx i?^"^/^, corresponding to the same total mass of the irradiated disk 
(assuming Rd = 373 AU), is also shown (dotted line) (see text). The disk and stellar parameters are the 
same than Figure ||. 
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Fig. 8. — Radial distribution of total mass of an irradiated (solid line) and a non irradiated (dashed line) 
disk model. The mass of the central star is Af* = 0.5 Mq. For the irradiated disk, the approximation given 
by equation (|3^) is also plotted (stars). The disk and stellar parameters are the same than Figure 
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Fig. 9. — Radial distribution of the viscous time for an irradiated (solid line) and a non-irradiated (dashed 
line) disk as a function of radial distance to the central star. The mean disk age estimated (Strom et al. 
1995) is ^ 10^ yrs. The disk and stellar parameters are the same than Figure ||. 
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Fig. 10. — Vertical distribution of the fraction of total energy flux transported by different mechanisms, for 
different annuli. The fluxes are: total intrinsic flux (solid line), radiative flux (dotted line), turbulent flux 
(dashed line) and convective flux (dot-dashed line). The position of the gas scale height evaluated at the 
midplane temperature H is shown with a vertical dashed line and the position of the photosphere is shown 
with a vertical dot-dashed line. Each panel shows at the upper corner the radius of the given annulus in 
units of stellar radius. The disk and stellar parameters are the same than Figure 
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Fig. 11. — Radial distribution of Toomrc's instability parameter Qt for an irradiated (solid line) and a 
non-irradiated (dashed line) disk model. The value Qt = 1 corresponds to the limit between stable and 
unstable regions (dotted line). The disk and stellar parameters are the same than Figure ||. 
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Fig. 12. — The contour of xjxi = 1, where x is the disk ionization fraction and Xi is the ionization fraction 
corresponding to a magnetic Reynolds number Rcm = 1- The region inside this contour, 0.2 < ii < 4 AU 
and z < 0.17, has Rcm < 1- We also plot the disk height Zoo (solid line), the height where a large fraction 

of the stellar radiation is absorbed Zg (dotted line), the photospheric height Zphot (dashed line) and the gas 
pressure scale height evaluated at the midplane H (dot-dashed line) 
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Fig. 13. — SED of the reference model pole-on, compared with AA Tau observed SED. The system's SED 
(star+disk) is shown with a solid line, the central star SED, taken from Bruzual & Chariot (1993), is plotted 
with a dotted line, and the disk SED is plotted with a dashed line. Also, the SED of a purely viscous disk 
with the same M, a and central star, is shown with a dot-dashed line. The observational points (squares) 

were taken from Adams et al. (1990), Bcckwith et al. (1990), Beckwith & Sargent (1991), Weaver & Jones 
(1992), Kenyon & Hartmann (1995), Dutrey et al. (1996). The error bars are related to the variability of 
the object for A < 12 /zm, and represent the reported errors for A > 12 /x. 
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